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Abstract: This paper is a theoretical analysis of mirror tilt in a Michelson interferometer 
and its effect on the radiant flux over the active area of a rectangular photodetector or 
image sensor pixel. It is relevant to sensor applications using homodyne interferometry 
where these opto-electronic devices are employed for partial fringe counting. Formulas are 
derived for radiant flux across the detector for variable location within the fringe pattern 
and with varying wave front angle. The results indicate that the flux is a damped sine 
function of the wave front angle, with a decay constant of the ratio of wavelength to 
detector width. The modulation amplitude of the dynamic fringe pattern reduces to zero at 
wave front angles that are an integer multiple of this ratio and the results show that the 
polarity of the radiant flux changes exclusively at these multiples. Varying tilt angle causes 
radiant flux oscillations under an envelope curve, the frequency of which is dependent on 
the location of the detector with the fringe pattern. It is also shown that a fringe count of 
zero can be obtained for specific photodetector locations and wave front angles where the 
combined effect of fringe contraction and fringe tilt can have equal and opposite effects. 
Fringe tilt as a result of a wave front angle of 0.05° can introduce a phase measurement 
difference of 16° between a photodetector/pixel located 20 mm and one located 100 mm 
from the optical origin. 
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1. Introduction 



The Michelson interferometer [1] shown in Figure 1 has been used extensively in the field of 
metrology, most famously in the Michelson-Morley experiment [2]. The measurement is based on 
detecting and measuring the number of complete and partial fringes resulting from translation of one of 
the reflecting mirrors and relating the result to the light-source wavelength. 

Figure 1. Michelson interferometer. 



Variable width area over 
which Radiant flux is derived 




Plane of Photodetector 



Mirror Ml (fixed) 



Beamsplitter (50:50 Lossless) 



^^^B Mirror M2 (movable, can be tilted) 

Discrete photodetectors and image sensors are commonly used to detect the sinusoidal pulsing 
fringe pattern. The radiant flux of the fringe pattern incident on the device active area is 
mathematically derived by integrating the irradiance over the circular aperture of the light source [3]; 
over a circular aperture of the interferogram [4-10]; over a square/rectangular aperture of the 
interferogram [8,11]. Effectively, the active area of the photodetector performs the same function on 
the irradiance giving an output proportional to the radiant flux. 

When using a well coUimated beam and plane flat mirrors that are not perfectly aligned, i.e., tilted, 
fringe lines of equal inclination, width and spacing are produced that contract as the tilt angle is 
increased and expand as the tilt angle is reduced, having a significant effect on the radiant flux over the 
active area. This change in radiant flux can be a source of measurement error [3-9] when unconsidered 
in applications [12-16]. 

As wave front angle increases, the modulation amplitude of the dynamic fringing reduces and at a 
specific tilt angle the modulation amplitude of the radiant flux becomes zero [3,5-11,17]. The 
modulation amplitude is found to decay as a cardinal sine function with a rate of decay proportional to 
the area of the photodetector and inversely proportional to wavelength. 

To overcome or minimise mirror tilt prevalent with flat plane mirrors, comer cube retro-reflectors [8,9], 
cat-eye reflectors [18] or alternate interferometers [11] are used. 

The behaviour of the radiant flux with varying wave front angle is also affected with varying 
photodetector distance from the central axis of the interferometer and varying its distance from the 
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optical model origin. Analysis of the radiant flux with varying wave front in conjunction with 
photodetector area, distance from beam centre, distance from origin and wavelength appears not 
to be covered in the literature although [17] has included distance of the photodetector from the 
optical model origin specifically to determine the maximum offset angle for a beam-tilting spatial 
modulation interferometer. 

Despite modulation amplitude vs. wave front angle being well understood, to the best of our 
knowledge, the following analysis is not covered in the literature, i.e., the behaviour of the radiant flux: 
for a rectangular aperture for varying wave front angle beyond the first modulation zero; over varying 
active area widths for varying distances from beam centre, e.g., row of pixels across an image sensor; 
for varying active widths and distances from the tilted mirror; with fringe contraction speed and mirror 
tilt; with decay constant. 

This paper addresses these issues, specifically: 

• Behaviour of the radiant flux for variable wave front angle as a function of photodetector 
width and position within the fringe pattern; 

• Behaviour of the radiant flux on two identical photodetectors adjacent each other; 

• Magnitude of the radiant flux at wave front angle(s) of equal radiant flux; 

• Polarity reversal of the radiant flux beyond specific wave front angles; 

• Behaviour of the radiant flux for variable wave front angle with variable distance of the 
photodetector from the tilting mirror; 

• Speed of transition of the fringe lines across the photodetector for variable wave front angle; 

• Speed of fringe line tilt across the photodetector for variable wave front angle; 

• Damping function constant of the radiant flux for variable wave front angle. 

The relevance of this theoretical analysis is to make evident how these other factors may have an 
undesirable effect on sensor applications using homodyne interferometry where photodetectors or 
image sensors are employed to sense small fractions of a fringe to achieve extremely high resolutions 
of measurement. It goes beyond the adverse effect of modulation amplitude reduction due to increasing 
wave front angle [3,5-11,17] and introduces what have been termed primary nodes. The results from 
the mathematical analysis describes how the radiant flux behaves when the five parameters; wave front 
angle, wave length, photodetector width and position (x, y) are varied independently and concurrently, 
and how this behaviour can introduce fringe counting errors. 

2. Mathematical Analysis 

The analysis is carried out based on a conventional Michelson Interferometer that is configured as 
shown in Figure 1 with the following configuration constraints: 

• Light source is a coUimated monochromatic beam; 

• Wave fronts over the area of the photodetector are approximated to be plane waves; 

• Flat plane mirrors are used to reflect the transmitted and reflected beams back to the beamsplitter; 

• Beamsplitter is lossless and is non-polarising and creates a transmitted and reflected beam of 
equal amplitude. 
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The 
follow 



two wave fronts are orientated as depicted in Figure 2 and all calculations are based on the 
ng constraints: 

y axis is taken to be normal to wave front 1 ; 

Origin of the Cartesian coordinate system is the point at which the centre of the incident 
beam is reflected by mirror M2; 
Mirror M2 tilts only about the z-axis; 
Mirror M2 translates only along the )^-axis; 
Plane of the photodetector remains orthogonal to the y axis; 

Shape of the active area of the photodetector is rectangular with variable side length s in the 
X direction and fixed side length z (set to unity) in the z direction; 
Fringe pattern irradiates the entire active area of the photodetector; 

Output of the photodetector is assumed to be a 1:1 linear function of the incident radiant flux; 
Distance to the photodetector from mirror M2 is variable. 



Figure 2. Wave fronts 1 and 2 with Mirror M2 tilted at angle 612 about the z-axis. 
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The mathematical analysis is divided into the following subsections, the outcome of which is 
studied further in the Results section: 

1. Derivation of the equation for radiant flux from irradiance of the fringe pattern; 

2. Identification of specific wave front angles 6„ with invariable radiant flux for two displaced 
photodetectors of equal size active areas; 

3. Determination of the magnitude of the radiant flux at specific wave front angles On', 

4. Determination of the linear equation defining the profile of the fringe pattern in the x-y plane; 

5. Determination of the speed of the fringe lines with variable wave front angle 6; 

6. Determination of the damping function of the radiant flux with variable wave front angle 6. 

Note: The angle 0 in this paper refers to the angle that Wave Front 2 makes with Wave Front 1. The 
tilt angle of Mirror M2 is therefore 6/2 relative to Mirror Ml. 

2.1. Derivation of the Equation for Radiant Flux 

The electric field of a plane wave is given by Equation (1) [19]: 

E(r, t) = Eoe^^'^'-'^*) (1) 

where E is the time (t) dependent electric field, r is the unit vector of the electric field in 3 dimensional 
space, i.e., r = xx + yy + zz and x, y and z are unit vectors along the x, y and z axes, Eq is the vector 
amplitude of the wave, k is the wave vector where k = ku, where u is the unit vector defining the 
direction of propagation [20] and |k| = k = 2n/A, X is the wavelength of the light source, k is the 
wave number and w is the angular frequency of the wave. 

Figure 2 depicts the linear optical equivalent of the Michelson interferometer with the virtual source 
wave front approaching mirrors Ml and M2 from the top of the figure. With reference to the origin, 
the reflected wave fronts 1 and 2 from respective mirrors have wave vectors k^ and k2 : 

ki = ky (2) 

k2 = /c sin 0 X + /c cos B y (3) 

Also depicted in Figure 2, the source wave front travels a distance AJ further to Ml creating an 
optical path difference (OPD) between the wave fronts and a phase lag of kllS.d relative to wave front 2. 
The sum of the electric fields of wave fronts 1 and 2 is therefore: 

Esum(r, t) = Epe'^'^i-'-'^^Ad-a.t) + Eoe'^*^2T-a.t) (4) 
EsumCx.y.z) = £'Qe-i«t(eife(y-2Ad) + gi(fc(y cos e+xsmS))^ 

The irradiance / of an electric field is given by Equation (6) and is the radiant flux of the electric 
field delivered per area to a given surface with units Wm~^, i.e., radiant flux density: 

^ ^ J ^sum ' ^sum (6) 

where nju is the refractive index of the medium, c is the speed of light in vacuum, cq is the vacuum 
permittivity, and ^sitm the complex conjugate of Esum- 
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/ = 2 (^^Bl^^ Eo^{l + cos(/c(y -2M-ycosO-x sin £»))) 

Equation (7) indicates tiiat the irradiance at a point (x, y) in the fringe pattern created by wave fronts 1 
and 2 is dependent on the values of x and y, wave number k, which is a function of wavelength, the 
angle 6 between the wave fronts and the optical path difference 2Ad. 

If Equation (7) is integrated along the x-axis between arbitrary points xi and X2 and then multiplied 
by side length z in the z-direction to create an area across the photodetector, the solution is the radiant 
flux incident on a rectangle of side lengths X2 - xi = s and z. As mirror M2 is only tilted about the 
z-axis, variable z does not need to be included in the integration as it behaves purely as a multiplier to 
the solution of the integration along the x-axis. Therefore: 

O g = z ■ J / dx = 2 ( ^^'^^^ ^ Eo^z j (l + cos(/c(y - 2Ad - y cos 0 - x sin d))) dx 

I 1 ^2 



\ksmd 



(sin(/c(x sin 0 + y cos G —y + 2Ad)) + kx sin d) 



The radiant flux given by Equation (8) is expressed in Watts (W) and is the total radiant 
power of the interference beam incident on the defined rectangular active area of the photodetector. 
Atd = 0,Oe = 0/0 which is indeterminate, therefore applying L'Hopital's rule to the integral solution of 
Equation (8) for ^ ^ 0 returns: 

(sm(k(x sm9 +y cos 9 —y + 2Ad)) + kx sin 9) 

lim — — —r^ ^ = x(l + cos(/c2Ad)) (9) 

0^0 ksm9 

Therefore, as ^ 0, the radiant flux derived in Equation (8) tends to: 

Oe^e^o) = 2 (^^^) Eo\ ■ s(l + cos(/c2Ad)) (10) 

The radiant flux in Equation (10) is a maximum when cos(k2Ad) =l,i.e., when 2Ad = n/X, where n/ 
is an integer equivalent to the number of fringe lines and 2Ad is the optical path difference. Whenever 
the OPD is an integer multiple of the wavelength, the two wave fronts in Figure 2 are in phase with 
one another resulting in maximum radiant flux, i.e. : 

^e(0-O,2Ad=n;^A) = 2 (^!^^^) Eq Z ■ 2s (11) 

To demonstrate the behaviour of the radiant flux over differing integral boundaries. Figure 3 shows 
the radiant flux curves for two sets of integral boundaries that are equal in length with assigned 
variables defined as follows that have been substituted in Equation (10);y = 0m, A = 680 x 10"^ m 
therefore k = 9,239,978, Ad = 0 m, integral width s = 0.001 m, red curve integral boundaries 
X2 = 0.0005 m, Xl = -0.0005 m, blue curve integral boundaries X2 = 0.001 m, xi = 0 m. 

It can be seen from the Figure 3 that there are node points at half the normalised radiant flux that are 
cyclic, which have been termed primary nodes, and this phenomenon is explored further below. What 
is also noticeable is the two curves converge as ^ ^ 0 as predicted in Equation (10). 
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Figure 3. Normalised radiant flux vs. varying 6 with y = Q, X = 680 nm, At/ = 0, integral 
width s = I mm, for integral boundaries X2 = 0.5 mm & xi = -0.5 mm and X2 - I mm & 
jci = 0 mm. 

Normalised radiant flux vs. varying 6 for overlapping equal width active areas 
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2.2. Identifying Specific Wave Front Angles On with Invariable 0efor Two Sets of Integral Boundaries 

To analyse the effect of mirror tilt angle on two separate rectangular areas of equal size and 
determine the node points observed in Figure 3, consider only the integral solution of Equation (8). If 
we define the two intervals along the .jc-axis with upper and lower limits xi, xj & x^, Xi, such that 
xi-x\=xa,-XT, = s and substitute in the Equation (8) we get after simplification Equations (12) and (13): 



O oc [sin(/c(x. sin 0 + y cos 0 - y + 2Ad)) + /cx, sin B 

K Sin u 

— sin(A:(xi sin 0 + y cos 6 — y + 2Ad)) — /cx^ sin O] 

« _1_ [3,„(,(., Sin « + y cos e - y + 2A.)) + fc., sin 0 
K Sin u 

— sin(A:(x3 sin 0 + y cos 6 — y + 2Ad)) — kx^ sin O] 



(12) 



(13) 



To solve for 0 and x, let Equation (12) = Equation (13), applying trigonometrical identities and 
simplifying yields: 

^ks sin 9\ \ //c((x2 + x^) sin 0 + 2y cos 0 — 2y + 4Ad)\ 



sm 



]■ cos I - 



— cos 



'/c((x4 + X3) sin 0 + 2y cos 0 - 2y + 4Ad)' 



(14) 



To obtain the node points that satisfy Equation (14) for the two intervals defined above, i.e., 
X2 ~ Xl = X4 — x^ = s, values need to be assigned to these boundary limits. For example, let Xl, X2 & X3, X4 
be the two intervals depicted along the plane of the photodetector in Figure 2 with values defined as 
Xl = -s; X2 = 0; X3 = -s/2; xa = s/2. 



Sensors 2013, 13 



11868 



Substituting these values in Equation (14) and simplifying yields: 

fks sin e\ \ fkU-s) sine + 2y cose -2y + AMY 
sin(^— j-|cos(^ 2 ^ 

//c(2y cos 0 - 2y + 4Ad)^ 



— cos 



(15) 



= 0 



The equality of Equation (15) is satisfied if: 

(a) sin ^^Ei^IL^j = q and/ or 

(b) COS (^^ j - COS J = 0 

Solving: 

(a) is satisfied when (ks sin ff)/2 = ripju, where Hp is an integer and is the number of what is 
termed a primary node (see Figure 3), resulting in: 

smen =^-oren =— (16) 

If small angles are considered, implying sin = 0^^, where the intersection of the two radiant 
flux curves at Q^^ is called a primary node and 6^^ is called the primary node angle 

(b) is satisfied when 

fkU-s) sin e + 2ycose -2y + 4M)\ fk(2y cos e-2y + 4M)\ 
cos j = cos j (17) 

There are two solutions that satisfy Equation (17): 

i. (s) sin ^ = 0 therefore sin Bo = 0 where Oo = 0. As only small angles are of concern, 0 = nnis 
irrelevant. Note that = 0 is co-incident with primary node np = 0 from Equation (16). 

ii. The cosines are identical for nsTc ± d/2, where ^ is a phase shift and ns is an integer related to 
secondary nodes, i.e.,: 

cos ^n^TT + =cos ^n^TT — (18) 

ns can be zero if the data is mirrored about ^ = 0 and = 1 if the data is mirrored about d = 7u. 

The occurrence of secondary nodes is unique and specific to the defined integral boundaries xi,X2& 
X3, X4, the values of k, y and A d. Consequently, the cosine expression in Equation (14) has to be solved 
accordingly with its own set of boundary conditions. The wave front angle at secondary node angles 
e^^ is incidental, unlike at primary nodes, which is cyclic and dependent only on k and s (sine 

expression in Equation (14)). 



2.3. Determination of the Magnitude of0e Wave Front Angles 9j 



p 



To determine the magnitude of the radiant flux of the fringe pattern at the primary nodes we solve 
Equations (12) and (13) independently for the intervals previously defined i.e., x\ = -s; X2 = 0; 
X3 = -5/2; X4 = s/2. Beginning with Equation (12): 
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sin(/c(y cos 9 — y + 2Ad)) — sin(/c(— 5 s\nd +y cos 0 — y + 2Ad)) + ks sin 0 

k^e ^'^^ 

For small angles sin 6 = 6, therefore Equation (19) becomes: 

sin(/c(y cos 0 - y + 2Ad)) - sin(/c(-s0 + y cos 0 - y + 2Ad)) + ksG 

At 6^, Equation (20) reduces to 0/0 , which is indeterminate. Therefore applying L'Hopital's rule: 

sin(/c(y cos 0 - y + 2Ad)) - sin(/c(-s0 + y cos 0 - y + 2A£i)) + ksd 

lim — V:; 

0-0 kd (21) 

= s(cos(/c2Ad) + 1) 
Repeating the above for Equation (13) yields the same result. 

Equation (21) shows that as ^ 0 the radiant flux becomes equal for the integral boundaries 
xi = -s; X2 = 0; X3 = -5/2; = s/2. Substituting the result of Equation (21) into Equation (8) gives 
Equation (22), which shows the magnitude of the radiant flux as ^ 0 is dependent on k and Ad and 
independent of y: 

^e(0-o) = 2 (!^^) Eo\ ■ s(cos (/c2Ad) + 1) (22) 

Note Equations (10) and (22) are equal, resulting in maximum radiant flux for translations where 
2Ad = nA: 

/ni{i6oC\ 2 

^e(0-O,2Ad=nA) = 2 ) Eq Z • 2s (23) 

To work out the value of radiant flux for all other values of 9^^ , i.e., ^1,2,3,..., substitute 
9^^ = UpA/s from Equation (16) into Equation (20), which renders: 

'^e(x„X2) °^ 5 (24) 

Repeating the above for Oe (xs, X4) again yields the same result, therefore the radiant flux at primary 
node angles 6^^ — 61,2,3,.. is: 

^<W) = 2(^)^0^-5 (25) 

Equation (25) also shows that the radiant flux at = 6^ 2,3,.. is half maximum (cf. Equation (23)), 
and in contrast to the radiant flux at 6^^ = given in Equation (22), '^6(0^23 ) independent of k 

and Ad. The reason for this is the width s of the active area is an integer multiple of the fringe line 
spacing at primary node angles = 6^ 2,3,...- 

When there is an exact multiple of fringe lines within the active area [18], the radiant flux across the 
active area is the mean of the maximum constructive and destructive interferences, i.e., 50%. This 
means that if the active area is moved in either direction along the x-axis, the radiant flux remains 
static at 0.5 normalised magnitude. If mirror M2 translates, there will be no change in the radiant flux 
despite the fringe pattern moving back and forth across the active area. 
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2.4. Effect of Distance x and y of Photodetector from Origin with Varying 6 

To deteraiine the effect of distance y of the photodetector from the origin with varying 6, the 
position and slope of the fringe lines needs to be determined. This is done by finding the instances of 
maximum value of irradiance in Equation (7), i.e., when: 

cos(/c(y — 2Ad — ycos0 — A:sin0)) = 1 (26) 

Equation (26) is true when: 

/c(y — 2 Ad — y cos 0 — X sin 0) = 2n^7r (27) 

Where n/is the n''^ fringe line. Solving for y gives: 

xsmd nfA + 2Ad 
(l-cos0)"^(l-cos0) ™ 

Equation (28) defines the profile of the fringe pattern in the x-y plane as illustrated in Figure 2, 
where sin 6/(1 - cos^) is the slope of the fringe lines and (n/A + 2Ad)/(l - cosff) is the y intercept. As 
only small angles are being considered it becomes indeterminate at ^ = 0, which stands to reason as the 
fringe pattern is uniform across the active area as well as the x-y plane and no fringe lines are present. 

By solving Equation (28) for n/, the number of fringes lines passing over a given point (x, y) can be 
calculated for 6 increasing or decreasing from zero to a given wave front angle: 

y(l — COS0) xsinO 2Ad 
nr J —-— (29) 

The y term in the above equation is positive for ^ 0 and is symmetrical in shape as a function of 6. 
For small angles sin 6 = 6, therefore the coefficient of x is a linear function of angle 6. The polarity of 
nf, is therefore dependent on the magnitude and sign of x, 6 and Ad. 

nf is the number of fringes counted at a point {x, y) as 6 is varied. Assume Ad = 0, there is a special 
case in Equation (29) when: 

y(l — cos 0) X sin 9 



(30) 



and the fringe count = 0 despite ^ > 0. In this case, fringe lines would have moved over point (x, y) 
in one direction as 6 is increased and then back again as 6 is increased further to the angle 6 that 
satisfies Equation (30). 

As discussed with Equation (28), the slope of the fringe lines is given by: 



y smQ (0 



(31) 



X (1 — cos 9) 

Where 612 is the angle of the normal of mirror M2 relative to the y-axis (Figure 2). Therefore, there 
is a set of points {x, y) coincident with the mirror normal that renders «/= 0. 

2.5. Deriving the Speed of Fringe Movement 



If Ad is considered static {Ad = 0), Equation (29) reduces to: 
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y(l — cosd) xsind 
n,=^— — (32) 

Taking the derivative of the x-terai delivers the speed of contraction/expansion of the fringe lines at 
the point (x, y): 

dUf X 

dsmd " ~X ^^^^ 

The speed of sideways deflection of the fringe lines at point {x, y) results from calculating the 
derivative of the y-term: 

dUf y sin 0 y 

:r^ = TT . =T^^"^ (34) 
dsmd /I Vl — sin^ d ^ 

The overall fringe movement speed is therefore: 

dUf X y 

;r^ = -T + Ttan0 (35) 
a sin 0 A A 

From Equation (35), the direction of fringe movement reverses if: 

X y 

-- + -tan0 = O (36) 

and therefore the angle of movement reversal Orev is: 

0™ = tan~^ - = sin~^ 



That is, when the perpendicular of the wave front from the origin points towards +x. The mirror tilt 
angle is therefore If y 0, drev ~^ ?r/2, i.e., the smaller y, the larger is drev- 

From Equation (29), the fringe count returns to zero if the numerical value of x- and y-terms cancel 
each other out, i.e.,: 

y — ycosd — xsin9 = 0 (38) 

Which has solution: 

2xy 

sin 9 = r 2 ^ 2^ (39) 
{x^ + y^) ^ ' 



This occurs at the angle On = o: 



2xy 

^n=o = sin-i (^^2+y2) (40) 



If X = y, then d„ = o = nil, if x > y, then 6n = o> ?r/2 and conversely, if y > x, then 6n = o< nil. 
The relationship between ^„ = o and drev, results in the identity of d^^Q = 29rev 
From Equations (37) and (40): 

1 . 2xy X 
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Resulting in: 



2xy 

sin 29fev = 2 sin Of-ev cos ^^ev — T^~, — 7T (42) 



and proving tlie identity of 0„=o = ^^rev being correct. 
2.(5. Determination of the Damping Function of the Radiant Flux Curve 

From Figure 3 it is evident that the normalised radiant flux curve oscillates about a level of 
0.5 radiant flux and decreases its wave amplitude about this level with increasing wave front angle 6. 
As this behaviour constitutes a damped function, the radiant flux decay function is derived subsequently. 

For a centred photodetector of width s and positioned at x = 0, with integral boundaries X2 = +s/2 
and xi = -s/2, y = 0 and A = 0, the radiant flux 'i>e(xi,x2) Equation (8) reduces to: 

„ . rks sin 6\ 
2sm(^— ) 

^^^^--^^^^ k^9 + ^ 

In Equation (43), the maximal radiant flux corresponds to 2s (cf. Equations (11) and (23): 
radiant flux = constant * z * 2^; z is considered unity as the z-direction is perpendicular to fringe lines). 
Normalising the radiant flux to 2s delivers the normalised radiant flux On: 

si n J ^ 1 

sk sin 6 2 

As the normalised radiant flux ranges from 0 to +1, it is converted to a range from -1 to +1 in order 
to compare it to a standard sine function: 

'ks sin 6 \ . cks sin 6 



ks sin 6^ 

<D„cx - ^ +1 (44) 



. //cs smy\ . //cs smy\ 

sm( — ^ 1 sm( — ^ 1 /'ks sin 9\ /ks9\ 

*n 2 — ^^—^ = — , - ^ = sine — « sine -— (45) 

" sk sin 9 ks sm 9 \ 2 ) \ 2 ) ^ ^ 

2 

for small 6. Considering that k = 2nlX and that a damped function (e.g., cardinal sine) corresponds to an 
undamped function times a decay function, the decay function Fd is identical to the reciprocal of the 
denominator {ks sin 0)12 in Equation (45). For small 6: 

2 1 XI 

'"^Tk-e^Ts-e 

Equation (46) is the constitutive equation of the radiant flux decay function at x = 0 and y = 0. Fd is 

2 A 

a reciprocal function of d, with a decay constant Cd of — or — . 

SK. TtS 

3. Results 



The radiant flux and the fringe count are influenced by a range of parameters, namely the wave 
front angle 6, the position x of the photodetector with respect to the centre line, the laser wave length X, 
the distance y between the mirror and the photodetector, and the side length s, i.e., the size of the 
photodetector, all of which are variable. The influence of these parameters is explained in a systematic 
way based on the equations derived in the Mathematical Analysis section. 
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3.1. Influence of 6 on the Radiant Flux 

The effect of 6 on the magnitude of the radiant flux also depends on the values of the other 
abovementioned variable parameters. In order to demonstrate this, the radiant flux is examined with 
four different conditions, in which some parameters are kept constant whereas others are variable. 

3.1.1. 6 = Variable, x = 0,y = 0, s = Variable 

Figure 4a shows the radiant flux curves of two different photodetector areas. If ^ = 0, the radiant 
flux is 100%. With increasing 6, the radiant flux decreases first, reaches the first primary node at a 
radiant flux magnitude of 50% and oscillates about the 50% level with decreasing radiant flux 
amplitude. The normalised radiant flux curve corresponds to a damped sine function with a damping 
function of 2/{sk6) or Xl{ns6). Multiplying a sine function of the form sin {nsOIX), where X/s represents 
the reciprocal value of the first node angle 6y^^ , delivers a damped sine wave, which, after adding 1 and 

dividing the sum by 2, results in the normalised radiant flux curve. In contrast to a standard sine wave, 
where the first minimum is at (3/2) 7i, i.e., at 1.5 6^^^, the non-linear decay rate causes the first 

minimum to be located at an angle of 1.4304 By^^, with a normalised radiant flux magnitude of 
0.3913832. This magnitude is independent of s and /I. 

Figure 4. Normalised radiant flux against wave front angle; (a) node points; (b) wave front 
angle at 99% radiant flux; = 10 |am, 100 |j,m, and 1 mm; X = 680 nm. 
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According to Equation (16), the smaller s, the larger is the wave front angle at the first node point. 
If 5 = 100 [im and 1 = 680 nm, the first and tenth primary nodes are at 9 = 0.3896° and 3.896°, 
respectively, and if 5 = 10 |um, the first primary node is at 3.896° (Figure 4a). The smaller is s, the 
slower the radiant flux decreases with increasing wave front angle. If * = 10 |a,m, 100 |a,m, and 1 mm, 
the angle at 99% radiant flux is located at 6* = 0.43°, 0.043°, and 0.0043° (Figure 4b), respectively. 
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According to Equation (16), s and X have opposite effects: reducing ^ by a factor of two results in the 
same angles of primary node points and 99% radiant flux as does a two-fold increase of X. Figure 5 
exemphfies this principle in a contour plot of equal angles of 99% radiant flux as a function of s and X. 

Figure 5. Contour plot of equal wave front angles of 99% radiant flux as a function of the 
side length 5 of a square photodetector area and laser wave length X. 




At this point it has to be mentioned that the same radiant flux curves are obtained if y > 0 and 
X = 3;(l/sin^ - 1/tan^) according to Equation (31). Equation (31) provides the solution for n/= 0, which 
deviates in x-direction if y > 0. This is explained in more detail below in the section dealing with 
influence on fringe count. 

3.1.2. 0 = Variable, x = Variable, y = 0, s = Variable 

Photodetector positions x ^ 0 changes the curves shown in Figure 4a insofar as the number of 50% 
radiant flux transitions is larger than the number of primary node points. The larger x, the more the 
radiant flux curve oscillates between the node points. In Figure 6, the radiant flux curve dX x = sl2 
intersects the 0.5 radiant flux level once between each pair of node points, the curve dX x = s does so 
twice, at .x: = 2^ four times and oX. x = \0s twenty times. The larger x, the higher is the density of the 
radiant flux curve filling up the area under the radiant flux curve at jc = 0 (Figure 6b), which acts like 
an envelope curve for radiant flux oscillations at larger x. This is insofar important to note as it shows 
that the modulation amplitude of the radiant flux across x (Figure 7), is unaffected by x. 
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Figure 6. Normalised radiant flux against wave front angle at ;c = 0, sll, s, 2s (a); and 
x = Q and 10* (b); s = 0.001, A = 680 nm; secondary nodes are marked with green dots in 
the first cycle up to the 2nd node point (a). 
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Figure 7. Normalised radiant flux across an jc-range of 4 mm at different wave front angles 
6 (in degrees); the dashed green lines indicate the x-position of the blue radiant flux curve 
shown in Figure 6a (5 = 1 mm, x = s/2, s. Is); the dashed purple line indicates the radiant 
flux level (intersections of green dashed lines and purple radiant flux curve) at the first 
triple secondary node (at 6 = 0.02597°); the 1st minimum refers to the radiant flux curve at 
x = 0 (Figure 6a). 
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Figure 6a shows secondary node points, i.e., intersections of the two curves at radiant flux 
magnitudes other than 50%. Independent of the position x, all curves intersect at the primary node 
points. The radiant flux at the primary node points is constant (50%), whereas the radiant flux at the 
secondary nodes is variable and a function of x. For example, at multiples of 0.02597° (Figure 6a), the 
three radiant flux curves of x = sl2, s and 2s, with s = 0.001, intersect; at the 3rd and 6th intersection, 
the secondary nodes are identical to the primary nodes (2nd and 4th). 

As the angle 0 increases, so does the number of fringe lines per unit x (Figure 7). At ^ = 0, the 
radiant flux is constant at 100%. After a slight increase in 6, the radiant flux oscillates between 100% 
and 0%, i.e., the maximal radiant flux is still very close to 100% (Figure 4b). Further increase in 6 
reduces the radiant flux amplitude, which fluctuates about 50% until the modulation amplitude 
converges to 0 at the 1st node point and remains constant at 50% radiant flux. Further increase in 6 
expands the modulation amplitude, however, the polarity of the radiant flux curve changes, i.e., peaks 
at X = 0 before the node point are converted to troughs after the node point. 

3.1.3. 9 = Variable, x = 0,y = Variable, s = Variable 

When introducing the distance y between the plane of the photodetector and the tilting mirror, the 
radiant flux curve can be entirely below or above the 50% radiant flux level, touching it only at the 
primary node points (Figure 8). The radiant flux curve is then superimposed by a further oscillation of 
a longer wave length. At the 6th primary node point of ^ = 0.23377°, the radiant flux curve does not 
cross the 50% radiant flux level; nevertheless, the polarity changes in the same way as shown in 
Figure 7. The distance y does not affect the primary node points according to Equation (16), whereas 
the secondary nodes are a function of y (as well as s and 6). 



Figure 8. Normalised radiant flux against wave front angle at two different x; X = 680 nm, 
y = 0.02 m. 
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Increasing y (Figure 9) has the same effect as increasing x (Figure 6): the radiant flux curve 
oscillates more frequently under the envelope of the radiant flux at y = 0 (Figure 9a). This does not 
affect the modulation amplitude of the radiant flux (Figure 9b), which remains the same across x at a 
specific angle 6, however, the centre fringe line is more deflected off-centre with 9, the larger is y 
(Figure 9b). 

Figure 9. Normalised radiant flux against wave front angle (a) and x (b); s = 0.00001, 
y = 0, 0.0002, 0.002, and 0.02 m; (b) shows the position of the centre fringe {i.e., fringe 
number 0) radiant flux at angle ^ = +1° and movement of the fringe pattern with increasing 
y (note that amplitude range and fringe density are independent of y). 




Wave front angle (deg) ' 
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3.1.4. 9 = Variable, x = Variable, y = Variable, s = Variable 



Figure 10 summarises the influence of 9, s, x, y and X on the normalised radiant flux. The difference 
between Figure lOa-d is that the radiant flux before the first primary node point decays slower the 
smaller s is. Figure lOd shows for small angles that the modulation amplitude remains constant across 
X and 9. 

Figure 10e,f shows with greater y, the more the centre fringe line deflects towards larger positive x 
and that fringe lines from the negative x-side cross over to the positive-side. 

The dotted lines in Figure 10a,b,e,f shows the primary nodes for s = \ mm, s = 0.5 mm and that the 
primary nodes are dependent on s and independent of y for constant A. 

Figure 10g,h shows with greater A for equivalent s (cf. Figure 10a,b) that the interval of primary 
node angles is greater. Also, the radiant flux decays slower before the 1st primary node for greater A. 

Figure 10g,h shows with greater A for equivalent s (cf. Figure 10a,b) that the radiant flux drops off 
slower before the 1st primary node. 

Figure 10 shows at ^ = 0°, the normalised radiant flux = 1 and is independent of x, s, y and A. 
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Figure 10. Normalised radiant flux against wave front angle ^ = 0° ... 0.2°) and distance 
from the beam centre {x = 0 mm ... 5 mm) at 3; = 20 mm, and X = 680 nm; (a): 5 = 1 mm & 
y = 20 mm, A = 680 nm; (b): s = 0.5 mm & 3^ = 20 mm, A = 680 nm; (c): ^ = 0.1 mm & 
y = 20 mm, A = 680 nm; (d): s = 0.01 mm & 3^ = 20 mm, A = 680 nm; (e): s = I mm & 
y = I m, A = 680 nm; (f): s = 0.5 mm & y = I m, A = 680 nm; (g): s = I mm & y = I m, 
A = 1,550 nm; (h): s = 0.5 mm & y = I m,A = 1,550 nm. 




(g) 



Sensors 2013, 13 



11879 



3.2. Influence ofx andy on the Fringe Count 

The number of fringe lines n/ passing over a point (x, y) within the fringe pattern is given by 
Equation (29), which is a function of the fringe lines tilting (y term), the fringe lines contracting/expanding 
(.x:-term) and moving mirror displacement (AJ). 

Assuming AJ = 0, as the moving mirror tilts from orthogonality, fringe lines are produced with a 
slope that is parallel to the normal of the mirror (Equation (31) and Figure 2). The 3^- term in Equation (29) 
is linked to the slope of the fringe lines, i.e., as the wave front angle increases then so does the fringe 
lines (in cross-section). With x and 6 static, the further the point is up the 3'-axis, the greater the number 
of fringe lines that tilt across the point. The angle 6 between the two wave fronts can be positive or 
negative with respect to the j-axis, however, the coefficient of the j-term is always positive and 
therefore only has an additive effect on nf. 

The x-term of Equation (29) is linked to contraction/expansion of the fringe lines with wave front 
angle. The focal point of the contraction/expansion is the normal of the mirror that is coincident with 
the axis of tilt. 

With AJ = 0, as ^ is increased a central fringe line is generated and aligns itself with this normal. 
Fringe lines develop to the left and right of this normal and contract toward it. For a given point (x, y) 
in the fringe pattern, as 6 is increased, more and more fringe lines will develop and cross over the point. 
As the central fringe line is essentially static, fringe lines to the far left and right move far quicker than 
those closer to the central fringe line. 

Figure 1 1 and Equation (33) confirm this behaviour showing that the speed of contraction/expansion 
(from the x-term of Equation (32)) of the fringe pattern is a constant. Therefore an active area located 
further from the centre of the beam will experience more fringe lines passing over it than one 
located closer to the centre when mirror M2 tilts. Parameters x and X have opposite effects on the 
contraction/expansion speed. 

Figure 11. The effect of variable y on the fringe count; (a) x = 0.001 m, y = 0.02 m, 
1 = 680 nm, Ad = 0; (b) x = 0.001 m, >; = 0.573 m, ^ = 680 nm, Ad = 0. 
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The speed of the fringe line tilt (from the j-term of Equation (32)) results from Equation (34). This 
fringe tilt speed is a tangent function of the wave front angle, independent of x, i.e., of the lateral 
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position of the photodetector, but dependent on y and X. The speed of fringe line tilt is initially smaller 
than the speed of contraction, as the former is zero if ^ = 0 (tan ^ = 0). 

Figure 12 shows the effect of increasing and then decreasing fringe numbers with progressive wave 
front angle. The fringe count on the positive x-side are acutely curved initially, the speed of the centre 
fringe line deflection lags behind the speed of contraction. Subsequently, the former speed term 
catches up and eventually overtakes the latter term. This results in the fringe lines initially moving over 
an off-centre photodetector in one direction and then moving over the same detector again but in 
opposite direction, thereby first increasing the fringe count and subsequently decreasing it. Figure 12 
also shows that the larger x, the larger is ^rev- 
Figure 12. Fringe counting as a function of x and 9; x = 0.001 m, y = 0.30 m, X = 680 nm. 
Ad = 0; the fringe count is shown on the right side, for x = 0, 0.001, and 0.0015 m; 
"0" = peak of centre fringe; positive and negative fringe numbers refer to negative and 
positive X, respectively, i.e., to left and right sides of the centre fringe line "0". 
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If Ad is dynamic and 0 is variable, then the fringe count is affected by transition of the fringe 
lines across the photodetector due to mirror translation in addition to fringe tilt and fringe 
contraction/expansion. The direction of transition of the fringe lines is dependent on whether 
translation of the mirror is in the positive or negative y direction. 

The effect of fringe tilt and y position of the photodetector can also be shown taking the 
photodetector parameters given in Figure lid where s = 0.01 mm, y = 20 mm, A = 680 nm and 
centering the photodetector on the y axis (i.e., x = 0 mm). From Equation (8) the normalised radiant 
flux is calculated to be 0.997 for a wave front angle d of 0.05°. Relocating the photodetector at 
y = 100 mm returns a normalized flux of 0.938, which is reduced from the first location as a result of 
fringe tilt. To work out what this change in radiant flux represents in terms of change in fringe position, 
substitute each flux value into the equation O = O ^ sin(/cx'), where = 1 is the maximum flux 
amplitude, k = Inldf, df = XI6 is the fringe width and x' represents the first and second position 
respectively of the fringe lines. Subtracting the two equations and solving for the change in fringe 
position gives 34.64 f^m. Fringe width is 779.2 \im, therefore the change in fringe position due to the 
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photodetector being located further away represents a phase difference of 16° and consequently 
represents a difference in n/ between the two positions with varying 6. 

4. Discussion 

The focus of this paper has been to establish the behaviour of the radiant flux of the interferogram 
over a photodetector of rectangular aperture that is variable; in size; in displacement across the 
interferogram; and, in axial distance from the source of interference, for variable angle between the 
two wave fronts and variable wavelength. 

The most apparent observation from the mathematical analysis in this study is that the radiant 
flux decays rapidly with increasing wave front angle with the recurrence of primary nodes where the 
radiant flux decays to 50% maximum and the modulation amplitude reduces to zero. This observation 
is also confirmed by [3,5-10,17] where radiant flux is calculated over a disc and also by [8,11] where 
the radiant flux is calculated over a square area. The modulation amplitude in this study and the 
literature is found to decay as a cardinal sine function. However, this study has gone further to 
determine that the radiant flux decays with a decay function that is a reciprocal function of the wave 
front angle with decay constant that is proportional to the wavelength and is inversely proportional to 
the photodetector active area width. 

What is so far not apparent from the literature, nor is it evident in the figures, is that between each 
primary node the polarity of the radiant flux reverses. This phenomenon has an adverse effect on phase 
measurement accuracy where the wave front angle has been allowed to increase beyond a primary 
node and the modulation amplitude is sufficiently large enough for fringes to be counted. 

In the literature [3,5-11,17], the boundary of the radiant flux calculation is centred on the 
interferogram, giving just a single radiant flux curve decaying as a cardinal sine function. Whereas, in 
this study, the radiant flux boundary is derived to be variable across the interferogram. This variant 
shows that as the boundary is moved off centre, the radiant flux oscillates increasingly about the 50% 
normalised maximum as an integer multiple of the distance from centre. Additionally, the amplitude of 
the oscillation is bounded by the centred radiant flux curve. 

A further finding from this study, which to the best of our knowledge, is not mentioned in the 
literature is how the radiant flux is affected when the fringe lines tilt and contract/expand with varying 
wave front angle. Having included as an integral parameter the axial distance of the photodetector 
from the interferometer, it is found for a centred photodetector that fringe tilt initially lags fringe 
contraction, but then fringe tilt becomes increasingly dominant on the radiant flux with increasing 
distance of the photodetector from the beamsplitter. 

Interferometry applications that use a plane flat mirror with translation stage, and discrete 
photodetectors [12-15] or position sensitive device [16] will suffer erroneous measurements if the 
wave front angle is not limited to give acceptable modulation amplitude. This can be done by choosing 
an appropriate photodetector aperture width that is an order of magnitude less than the fringe line 
spacing. Increasing the wavelength improves modulation ampUtude for equivalent photodetector 
aperture widths. 

If the wave front angle varies beyond primary node angles the radiant flux reverses polarity 
introducing an additional error with a k radian phase change in photodetector output. 



Sensors 2013, 13 



11882 



Alignment of the photodetector with the centre of the interferogram reduces the susceptibility of the 
fringe lines crossing over the photodetector as they contract as a linear function of distance from centre 
with increasing wave front angle. 

Finally, the tilt angle of the fringe lines increases with increasing wave front angle and they cross 
over the central axis of the interference beam with increasing distance from the beamsplitter. Therefore 
to limit the error in measurement that this produces, the photodetector should be located as close as 
possible to the beamsplitter. 

5. Conclusions 

The radiant flux O across the active area of a photodetector is a damped sine function (i.e., a 
cardinal sine function) of the wave front angle 6, with a reciprocal decay of the ratio of wavelength A to 
detector width s. The larger s and the smaller A, the faster is the decay of the radiant flux. 

If the radiant flux magnitude at wavefront angle ^ = 0 is normalised to 100%, then the radiant flux 
magnitude at the primary node points is 50% and at the first minimum is 39.14%. The polarity of the 
radiant flux changes exclusively at every primary node point. The radiant flux magnitude at specific 
6 is independent of any parameter if the centre of the fringe beam coincides with the centre of 
the photodetector. 

The larger the distance x of the photodetector from the centre of the fringe beam and/or the distance 
y between mirror and photodetector, the more the radiant flux oscillates under the envelope radiant 
flux curve generated if x and y = 0. These two parameters do not affect the modulation amplitude of 
the radiant flux. 

The movement of fringe lines with increasing (9 is a combined effect of fringe contraction 
(x-dependent; the faster the more the detector is off centre) and fringe tilt (j-dependent; the faster the 
larger y). Fringe contraction and tilt movement can have opposite effects, with fringe contraction 
lagging behind fringe tilt, such that the fringe count first increases, then decreases and then returns 
to zero. 

Consequently, significant fringe count errors occur if the photodetector is operated near or beyond 
primary nodes where radiant flux modulation reduces to zero and then changes polarity, or if the x and 
y distances of the photodetector are large therefore increasing the fringe contraction and fringe 
tilt influence. 
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Appendix 

This appendix provides a detailed derivation of the equations in the paper "Theoretical Analysis of 
Interferometer Wave Front Tilt and Fringe Radiant Flux on a Rectangular Photodetector" by authors 
R.M. Smith and F.K. Fuss. 

The equation numbers in this appendix follows the equation numbering in the parent paper, 
however, where additional equations are included they are designated an alphanumeric reference. 

Al. Derivation of the Equation for Radiant Flux 

The electric field of a plane wave is given by Equation (Al) 

E(r,t) = Eoe'^"-'-'^^^ (Al) 

where E is the time (t) dependent electric field, r is the unit vector of the electric field in 3 
dimensional space, i.e., r = xx + yy + zz and x, y and z are unit vectors along the x, y and z axes, Eq 
is the vector amplitude of the wave, k is the wave vector where k = ku, where u is the unit vector 
defining the direction of propagation and |k| = /c = — , A is the wavelength of the light source, k is the 

wave number and o) is the angular frequency of the wave. 

Figure 2 depicts the linear optical equivalent of the Michelson interferometer with the virtual source 
wave front approaching mirrors Ml and M2 from the top of the figure. With reference to the origin, 
the reflected wave fronts 1 and 2 from respective mirrors have wave vectors kj and k2 . 

fei = ky (A2) 

k2 = k sin dx + k cos 6 y (A3) 

Also depicted in Figure 2, the source wave front travels a distance Ad further to Ml before it is 
reflected. Therefore, the optical path difference (OPD) between wave fronts 1 and 2 is 2Ad resulting in 
wave front 1 having a phase lag of /c2Ad relative to wave front 2. 

The sum of the electric fields of wave fronts 1 and 2 is therefore 

_ £'^g-ia)t^gifc(y-2Ad) gt(fe(ycos0+xsin0))^ K^-^) 



The irradiance / of an electric field is given by Equation (A6) and is the radiant flux of the electric 

ith units 



field delivered per area to a given surface with units Wm ^. i.e., radiant flux density 



j_fnRifoC\ 

\ 2 ) ^""^ ^sum (Ad) 

Where n^; is the refractive index of the medium, c is the speed of light in vacuum, Eq is the vacuum 
permittivity, and ^sMm the complex conjugate of fi'sMTn- 

J _ ^^:^lf2^^^Q-io)t^Qik(y-2^d') ^. ^ikiy cos 6 +x sin 0)^ 

■ EqB^^^ {e~^^^^~^^^^ + g~^'^(>''^°s0+xsin0)'^ (A7) 

= Eq^ (2 + 2 cos(/c(y - 2Ad - y cos 0 - x sin 0))) 
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= 2 £"0^(1 + cos(/c(y - 2Ad - y cos 0 - x sin £»))) 

Equation (A7) indicates tiiat tiie irradiance at a point (x, y) in the fringe pattern created by wave 
fronts 1 and 2 is dependent on the values of x and y, wave number k, which is a fiinction of wavelength, 
the angle 9 between the wave fronts and the optical path difference 2Ad. 

If Equation (A7) is integrated along the x-axis between arbitrary points x^and X2 and then multiplied 
by side length z in the z-direction to create an area across the photodetector, the solution is the radiant 
flux incident on a rectangle of side lengths X2 — x^ = s and z. As mirror M2 is only tilted about the 
z-axis, variable z does not need to be included in the integration as it behaves purely as a multiplier to 
the solution of the integration along the :v:-axis. Therefore: 

Og = ^ ■ 1 1 dx = 2 (~^y^) z I (1 + cos(/c(y — 2Ad — y cos 6 — xsin 6))) dx 

' X. (A8) 



— - — J £"0 z ■ — — - (sin(/c(x sin 0 + y cos 6 —y + 2Ad)J + kx sin 6) 



The radiant flux Og given by Equation (A8) is expressed in Watts (W) and is the total radiant power 
of the interference beam incident on the defined rectangular active area of the photodetector. At 9 = 0, 
Oe = ^ which is indeterminate, therefore applying L'Hopital's rule to the integral solution of 

Equation (A8) for 0 ^ 0 returns: 

(sin(fc(A: sin 0 + y cos 9 —y + 2Ad)) + kx sin 9) 
ksm9 

(cos(/c(x s\n9 + y cos 9 —y + 2Ad))) ■ k(x cos 9 — ysin9) + kx cos 9 

~ e^o kcos9 (A9) 

(cos(/c2Ad)) • kx + kx 

= lim ; 

0^0 k 

= x(l + cos(/c2Ad)) 

Therefore, as 0 -> 0, the radiant flux derived in Equation (AS) tends to 

*e(0-o) = 2 (!^^) ■ |x(l + cos(/c2Ad))|^^ 
= 2 Eo z ■ s(1 + cos(/c2Ad)) 

It is worth noting that at 9 = 0, the irradiance in Equation (A7) reduces to 



(AlO) 



ke=o) = 2 Eo\l + cos(/c2Ad)) (AlOa) 



and differs with Equation (AlO) by the area z ■ s. 

Equation (AlO) indicates that as 0 ^ 0 the magnitude of the radiant flux is dependent on the area 
z • s,k and Ad but is independent of y. 
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The irradiance and radiant flux in Equations (AlO) and (AlOa) are maximum when cos(/c2Ad) = 1, 
i.e., when 2 Ad = UfX, where Uf is an integer equivalent to the number of fringe lines and 2 Ad is the 

optical path difference. Whenever the OPD is an integer multiple of the wavelength, the two wave 
fronts in Figure 2 are in phase with one another resulting in maximum radiant flux, i.e.. 



^e(6^0.2^d=nfX) = 2 ) Eq Z ■ 2s (All) 



2 

To demonstrate the behaviour of the radiant flux over differing integral boundaries, Figure 3 shows 
the radiant flux curves for two sets of integral boundaries that are equal in length with assigned 
variables defined as follows that have been substituted in Equation (AlO); y = 0 m, 1 = 680 x 10"^ m 
therefore k = 9239978, Ad = Om, integral width s = 0.001m, red curve integral boundaries X2= 0.0005 m, 
Xi= —0.0005 m, blue curve integral boundaries X2 = 0.001 m, Xi= 0 m. 

It can be seen from the Figure 3 what appears to be node points at half the normalised radiant flux 
that are cyclic, which have been termed primary nodes, and this phenomenon is explored further below. 
What is also noticeable is the two curves converge as 0 ^ 0 as predicted in Equation (AlO). 

A2. Identifying Specific Wave Front Angles 0^ With Invariable Ogfor Two Sets of Integral Boundaries 

To analyse the effect of mirror tilt angle on two separate rectangular areas of equal size and 
determine the node points observed in Figure Al, consider only the integral solution of Equation (A8). 
If we define the two intervals along the x-axis with upper and lower limits x-^^, Xi & ^^3, X4 
such that x^ — Xx= X/^— X'^= s and substitute in the Equation (A8) we get after simplification 
Equations (A12) and (A13): 

ksxxvd^ ^ ^ ^ ^ ^ ^ (A12) 

— sin(/c(xi sin 0 + y cos 0 — y + 2Ad)) — kx-^ sin Q\ 

1 

oc . [sin(/c(x4 sin 0 + y cos 0 - y + 2Ad)) + kx/s, sin B (A13) 
fC sin \j 

— s\xv{k{x-i sin 0 + y cos 0 — y + 2Ad)) — kx-^ sin Q\ 
To solve for Q and x, let Equation (A12) = Equation (A 13) 

sin(/c(x2 sin 0 + y cos 0 — y + 2Ad)) + kx^ sin B 

— sin(A:(xi sin 0 + y cos 0 — y + 2Ad)) — kx^ sin Q 

- sin(/c(x4 sin 0 + y cos 0 - y + 2Ad)) + /CX4 sin B (A13a) 

— sin(/c(x3 sin 0 + y cos 0 — y + 2Ad)) — /CX3 sin B 

sin(/c(x2 sin 0 + y cos B — y + 2Ad)) — sin(/c(xi sinB +y cos B — y + 2Ad)) 

= sin(/c(x4 sinB +y cos B — y + 2Ad)) — sin(/c(x3 sinB +y cos B —y + 2Ad)) (^l^^*) 

Applying the identity below and simplifying Equation (A13b) yields Equation (A 14). 
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2 sin 



/A -tl\ /A+ tl\ 
sin A — sin fi = 2 sin I — - — I cos I — - — I 

//c((x2 sin 0 + y cos 9 — y + 2Ad) — (x^ sind +y cos 0 — y + 2Ad))^ 



(A13c) 



cos 



^/c((x2 sin 0 + y cos 6 —y + 2Ad) + (xi sin 0 + y cos 9 —y + 2Ad))' 



— 2 sin 



^A:((x4 sin 0 + y cos 9 —y + 2Ad) — (X3 sin 0 + y cos 9 —y + 2Ad))^ 



(A13d) 



cos 



^/c((x4 sin 0 + y cos 8 —y + 2Ad) + (X3 sin 0 + y cos 9 —y + 2Ad))^ 



sin 



ks sin 0\ / /c((x2 + Xi) sin 0 + 2y cos 0 — 2y + 4Ad)^ 



^/cs sin y j 



cos 



^/c5 sin 0\ [k(^(.X4 + X3) sin 0 + 2y cos 9 — 2y + 4Ad)^ 



//C5Sint/\ , , ^ ^ ^ , 
- j 2 ' = ° 



sin 



^/cs sin 0^ 



cos 



^/c((x2 + Xi) sin 0 + 2y cos 0 - 2y + 4Ad)^ 



— cos 



^/c((x4 + X3) sin 0 + 2y cos 0 - 2y + 4Ad)^ 



= 0 



(A13e) 



(A14) 



To obtain the node points that satisfy Equation (A14) for the two intervals defined above, i.e., 
X2 — Xi = X4 — X3 = s, values need to be assigned to these boundary limits. For example, let x^, 
X2 & X3, X4 be the two intervals depicted along the plane of the photodetector in Figure Al with 
values defined as = —s; = 0; X3 = — s/2; X4 = s/2. 

Substituting these values in Equation (A14) and simplifying yields: 



sin 



ks sin 9\ fk(i-s) sin 0 + 2y cos 0 - 2y + 4Ad)^ 



^/cs sin y j 



cos 



ks sin 9\ //c((0) sin 0 + 2y cos 0 - 2y + 4Ad)' 



/KS sin u\ 



cos 



= 0 



(A14a) 



/ks sin 9\ 



cos 



'/c((-s) sin 0 + 2y cos 0 - 2y + 4Ad)' 



— cos 



^/c(2ycos0-2y + 4Ad)' 



(A15) 



= 0 



The equality of Equation (A 15) is satisfied if 

1. sin ( ^^ ^ j = 0 and/ or 

/fe((-s)sin0+2ycos0-2y+4Ad)\ / 

2. COS ( -^^ -j - COS ( 



fc(2 y cose-2y+4Ad) j _ ^ 
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Solving: 

(a) is satisfied when ^ = npTr, wliere Up is an integer related to what are termed primary 

nodes (see Figure 3), resulting in 

2Ti„7r 



sin d-n = ^^^^ or = — 
"p fcs "p s 



(A16) 



if small angles are considered, implying sin^^^ = 6^^, where the intersection of the two 
radiant flux curves at 9^^ is called a primary node and 9^^ is called the primary node angle. 

(b) is satisfied when 

//c((-s) sin e + 2ycose -2y + 4Ad)\ //c(2y cos 9 -2y + 4Ad)\ 
cos ^ j = cos j (A17) 

There are two solutions that satisfy Equation (A17): 

(i) . (— s) sin 9 = 0, therefore sin 9^ = 0 where 9q = 0. As only small angles are of concern, 

9 = nn is irrelevant. Note that 0o — 0 is co-incident with primary node = 0 from 
Equation (A16) 

(ii) . The cosines are identical for ngTi + -, where 5 is a phase shift and rig is an integer related 



to secondary nodes, i.e., 

cos {usn + 2) ~ i^s'^ ~ 2) 



(A18) 



rig can be zero if the data is mirrored about 0 = 0 and rij = 1 if the data is mirrored about 
9 = n. 

From Equations (A17) and (A18) there are 2 identities: 



cos/c— ssin0+2ycos0— 2y+4Ad2=cosns7r— 52 



cos 



'k{2ycos9 - 2y -I- 4Ad)' 



= cos {usTt + ^ 



As the left- and right-hand terms refer to the same angle 9 

ki—s sin 9 + 2y cos 9 — 2y + 4Ad)=2nsn — 5 
k(2y cos 9-2y + 4Ad)=2ns7r -I- S 

k(-s sin 9 + 2ycos9 -2y + 4M)+5=2ngn 

k{2y cos 9 -2y + 4Ad) - S=2nsn 

Therefore: 

k{-s sin 9 + 2ycos9 -2y + 4Ad)-l-5=/c(2y cos 9-2y + 4Ad) - S 
and 

25 = k(2y cos 9-2y + 4Ad) - k(-s sin 9 + 2ycos9 -2y + 4Ad) 
= k2y cos 9 - k2y + k4M + ks sin 9 - k2y cos 9 + k2y - k4M 
= ks sin 9 



(A18a) 

(AlSb) 

(A18c) 

(AlSd) 

(A18e) 
(AlSf) 

(A18g) 
(AlSh) 
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The IngTi term in Equations (A18e) and (A18f) lies exactly between 
'/c((-s) sin e + lycos e -2y + 4Ad)\ 



cos 



(A18i) 



and 



cos 



^{lycosB -2y + A^dy 



(AlSj) 



as revealed from Equation (A18h) 
Therefore: 



lUsU = k(2y cos 9 - 2y + 4Ad) - 



ks sin 0 



This eliminates the unknown term 5 and provides the function for secondary nodes. 
Equation (A 18k) renders: 

ZUsTT = k [ly cos B -2y + 4M - 



s sin 0 



2 

( ssin0\ 
UgU = k\v cos 6 —y + 2Ad — 1 



ky cos 9 — ky + k2Ad — 



4 

ks sin 6 



(A 18k) 
Simplifying 

(A181) 
(A 18m) 



Rearranging and grouping terms yields: 

ks sin 6 



= —n^n + ky cos 9 — ky + k2Ad 

n^TT 4ycos0 4y 8Ad 

sin 6 = 1 1 

ks s s s 



sin' 



^ = l^J ™^^ + ItJ +(— J -n^'^^'^^n^^ k^ 



32y2 



cos 6 + 



64yAd 



cosd 



64yAd 



sin' 



SZnjTry 64ns7rAd 



32y2 64yAd 64yAd 
—cos 9 H — cos 8 1 



/4n,7r\' /4y\' 2 a ^ /'^yx' /8Ad\' 32n,7ry 



32ns7ry 64ns7rAd 32y^ 



cos 0 



64yAd 64yAd 
+ „ cos e ^ 1 + 1 - sin^ 0 = 0 



64yAd 64yAd 
+ , cos 6 ^ 1 + cos^ 0 = 0 



+ 



32ns7ry 64ns7rAd 32y^ 



cos 6 



+ 1 1 cos2 e + 



( 64yAd 32y^ 32ns7ry^ ^ ^ 32ns7ry 64ns7rAd /4n7r\^ /"^^A^ 



+ 



/8Ad\^ 64yAd ^ _ 



ks- fes^ 



\ s / 



(A18n) 
(AI80) 

(A18p) 
(A18q) 



(A18r) 
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Let 

2 



^-{-f) +1 (A18s) 

64yAd 32y^ ^lUgiiy 
B = -J2 J2 ^ (A18t) 

(^n^W /4y\^ /8Ad\^ 32n„ Try 64n„ TrAd 64yAd 



then 

i4 cos^ e + Bcose + C = Q 



and 



(A18v) 



-fi + V52 - 4^C 

cos e = C0Si,2 = ^ (AlSw) 

By way of demonstrating the presence of secondary nodes for the case with integral boundaries 
defined above as = — s; X2 = 0; X3 = — s/2; X4 = s/2, firstly, assume Ad = 0 m. 
Equations (AlBs), (A18t) and (AlSu) now reduce to: 



(A18x) 

32y^ SZn^Try 



^2 1' c2 



s2 ks^ (AlSy) 
2 ./I,.. 2 



32ns7ry ^ 



ks^ (A18z) 

Then, the variables in Equations (A18x)-(A18z) are assigned values, for example those defined 
in Equation (A18A). 

3; = 0.02 m 

A = 680-10"'' m therefore k = 9239978 (A18A) 
integral width 5 = 0.001m 
Finally, Equation (A18w) is solved for different rig (positive and negative) to find where the 
radiant flux curves intersect at the secondary node points. 

Secondary nodes for the above defined conditions are returned for 

""'^ 2A 
when Ug < 0 and > 0 and when 0 < < 4 for < 0. 
Also 



_ -B - ^B'^ - 4AC 

when < 4 for 2 < 0 



cos 0„^2 = ^ (A18C) 
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Consider the solution defined by Equation (A18B) and let = —1. The secondary node point 
is calculated by this solution, i.e., 0_i = 0.07409° and confirmed by substituting the above 
conditions into the integral part of Equation (A8). 

Mirror M2 in Figure Al can tilt left or right of the y-axis, therefore 9^^ can be positive or negative, 
therefore returning secondary nodes of both polarities for in Equations (A18B) and (A18C). 

From the above derivation of the secondary nodes, a specific solution was obtained for the 
defined intervals specified in as = —s;x2 = 0; :x:3 = — s/2;x4 = s/2 and the variable values 
assigned in Equation (A18A). The occurrence of secondary nodes is unique and specific to the 
defined integral boundaries x^, Xi & x^^, X4, the values of k, y and Ad. Consequently, the cosine 
expression in Equation (A 14) has to be solved accordingly with its own set of boundary conditions, 
such as Equation (A18A). The wave front angle at secondary node angles is incidental, unlike 
at primary nodes, which is cyclic and dependent only on k and s (sine expression in Equation (A14)). 
The next section derives the magnitude of the radiant flux at the primary nodes. 

A3. Determination of the Magnitude of0g at Wave Front Angles 0„ 

To determine the magnitude of the radiant flux of the fringe pattern at the primary nodes we solve 
Equations (A12) and (A13) independently for the intervals defined as x^ = — s; X2 = 0; X3 = — s/2; 
X4 = s/2. Beginning with Equation (A12): 

1 

*e(xi,x2) ^ k^sinO [^'"(''■^^2 sin 0 + y cos G —y + 2Ad)) + kx2 sin 9 

— sin(/c(xi sin 0 + y cos 6 —y + 2Ad)) — kx^ sin 9] 

1 

oc ^^.^^ [sin(fe(0 sin 0 + y cos 0 - y + 2Ad)) + feO sin 9 (A19) 

— sin(fe(— s sin 0 + y cos 9 — y + 2Ad)) + ks sin 9)] 
sin(/c(y cos 9 — y + 2Ad)) — sin(fe(— s sin 0 + y cos 9 —y + 2Ad)) + ks sin 9 



oc 



ksm9 



For small angles sin 6 = 6, therefore Equation (A 19) becomes 

sin(fc(y cos 0 - y + 2Ad)) - sm{k{,-s9 + y cos 0 - y + 2Ad)) + ks9 
*«(^i'^2) ^ Y9 (A20) 

At Bq, Equation (A20) reduces to ^, which is indeterminate. Therefore applying L'Hopital's rule 

sin(fc(y cos 0 - y + 2Ad)) - sin(fc(-se + y cos 0 - y + 2Ad)) + kse 
lim — 

e^o kO 

cos(/c(y cos 0 — y + 2Ad)) ■ (— fey sin 0) — cos(fe(— s0 + y cos 0 — y + 2Ad)) ■ (fe(— s — y : 

= lim — ; — 

fl^o fc 

cos(fe(y - y + 2Ad)) ■ (-fcyO) - cos(fc(-sO + y - y + 2Ad)) ■ (fc(-s - yO)) + fes 
"e™ fc (A21) 

- cos(fc(2Ad)) ■ (fc(-s)) + fcs 
~ 61-.0 fc 

fcs cos(fc(2Ad)) + fcs 
~ 0^0 k 
= s(cos(fe2Ad) + 1) 
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Repeating the above for Equation (A 13) 

^e{x^.x^) « fcie [^'" ^ (I sin 0 + y cos 0 - y + 2Ad)^ + /c| sin d 

— sin ^/c ^— I sin 0 + y cos 6 — y + ZAdj^ + /c| sin 0) j (A21a) 

sin ^fc(|sin0 + y cos0 - y + 2Adfj + ks sin 6 — sin (k(-^smG + y cos0 - y + 2Ad)^ 

fesin0 

For small angles sin 0 = 9, therefore Equation (A21a) becomes 

sin (/c(|0 + y cos 0 - y + 2Ad)) + ksO - sin (k{-^d + y cos 0 - y + 2Ad)) 

^e(x3,x.) « ^ ^ ^ (A21b) 

At Equation (A20) reduces to ^, which again is indeterminate. Applying L'Hopital's rule 

sin (k{^9 + y cos 0 - y + 2Ad)) + ks9 - sin (k{-^9 + y cos 0 - y + 2Ad)) 

Um — 

8^0 k9 

cos (fe(f 0 + y COS e - y + 2Ad)) ■ fc(f - y sin e) + fcs - cos (fc(-f 0 + y cos 6 - y + 2Ad)) ■ fc(-| - y sin 0) 
~ k 

cos(/c2Ad) ■ + fcA: - cos(fc2Ad) ■ /c(-|) ^^^^^^ 
~ 0^0 k 

ks cos(/c2Ad) + ks 

= lim ^- — = s(cos (/c2Ad) + 1) 

Equations (A21) and (A21c) show that as 0 ^ 0 the radiant flux becomes equal for the integral 
boundaries defined as = —s ; X2 = 0 ; x-^ = — s/2 ; 3:4 = s/2 . Substituting the result of 
Equation (A21) or Equation (A21c) into Equation (A8) gives Equation (A22), which shows the 
magnitude of the radiant flux as 0 ^ 0 is dependent on k and Ad and independent of y. 

fn^ieQC\ 2 

^eie^o) = 2 ^ }Eo z- s{cos (klM) + 1) (A22) 

Note Equations (AlO) and (A22) are equal, resulting in maximum radiant flux for translations 
where 2Ad = nX. 

/ni{i6oC\ 2 

^e(e^0.2Ad=nX ) = 2 — j Eq z • 2s (^23) 

n X 

To work out the value of radiant flux for all other values of 6^^, i.e., 01,2,3,.., substitute = 

from Equation (A 16) into Equation (A20) and Equation (A21b) independently. Beginning with 
Equation (A20) and using the identity in Equation (A13b): 

sin my cos - y + 2Ad M - sin m -s^ + y cos - y + 2Ad M + /cs^ 

oc (A23a) 

, nA 
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sin cos ~ y + 2Ad^^ — sin |^/c ^— n^A + y cos ~ y + 2Ad^^ + /cnA 



oc 



kUpA 



2 sin 



/c (y cos - y + 2Ad j - k (^-UpA + y cos -y + 2Ad^ 

2 

k (^y cos (^^^ - y + 2Ad^ + k (-ripA + y cos - y + 2Ad^ 



oc 



kripX 
s 



2 sin 



^^^j ■ cos [k ^y cos - y + 2Ad - + /cripA 



oc 



kUpA 



2 sin^ripTrj ■ cos I /c(y cos I I — y + 2Ad ) I + kripX 



oc s 



Repeating the above for Equation (A21b) 



Sin m f ^ + y cos - y + 2Ad M + fcs ^ - sin m - 1^ + y cos - y + 2Ad) 
oc z 



sin 



oc 



oc ■ 



k + y cos - y + 2Ad^ j + fcripA - sin [k{-^ + y cos (^^j - y + 2Ad^ 



knX 
s 



kUpA 



ky-^ + y cos (-^) - y + 2AdJ - fc — 2- + y cos (-^^ - y + 2Adj 



2 sin 



cos 



(A24) 



/c(^ + ycos(^)-y + 2Ad) + /c(-^ + ycos(^)-y + 2Ad) 



oc 



s 



oc 



kUpA 



2 sin ■ cos ^y cos ~ y + 2Ad j + fenp/lj 

2 sin(np7r) ■ cos ^y cos ~y + 2Ad j + fcripA 



oc s 
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Equations (A23a) and (A24) give the same result for the integral of the active areas at 9^^ = ^1,2,3,.. 
and therefore the radiant flux at these wave front angles is given by Equation (A25). 



x.,3,...) = 2(^)^0^-5 (A25) 



e(ei,2,3,...) - ^ V 2 

Equation (A25) also shows that the radiant flux at 6-^^ — 0i,2,3,.. is half maximum 
(cf. Equation (A23)), and in contrast to the radiant flux at 9^^ = 9q given in Equation (A22), 
^e(ei2 3 ) independent of k and Ad. The reason for this is the width s of the active area is an integer 
multiple of the fringe line spacing at primary node angles 9^^ = 0i,2,3,..- 

When there is an exact multiple of fringe lines within the active area, the radiant flux across the 
active area is the mean of the maximum constructive and destructive interferences, i.e., 50 %. This 
means that if the active area is moved in either direction along the x-axis, the radiant flux remains 
static at 0.5 normalised magnitude. If mirror M2 translates, there will be no change in the radiant flux 
despite the fringe pattern moving back and forth across the active area. 

A4. Effect of Distance x andy of Photodetector from Origin with Varying 9 

To determine the effect of distance y of the photodetector from the origin with varying 0, the 
position and slope of the fringe lines needs to be determined. This is done by finding the instances of 
maximum value of irradiance in Equation (A7), i.e., when 

cos(/c(y - 2Ad - y cos 0 - x sin 0)) = 1 ^^^^^ 
Equation (A26) is true when 

/c(y-2Ad-ycos0-xsin0) = 2n/7r ^^^'^^ 

where Uf is the n*'* fringe line. Solving for y gives: 

y (1 — cos 9) = xs\n9 + UfX + 2Ad 

(A27a) 

X sin 9 UfX + 2Ad 
^ " (1 - cos 0) (1-COS0) (A28) 

Equation (A28) is a linear equation and defines the profile of the fringe pattern in the x-y plane as 

sine nfX+2hd 

illustrated in Figure Al where ; — is the slope of the fringe lines and — is the y intercept. As 

(1— COS 0) (1— COS Q) 

only small angles are being considered it becomes indeterminate at 0 = 0, which stands to reason as 
the fringe pattern is uniform across the active area as well as the x-y plane and no fringe lines are present. 
By solving Equation (A28) for Uf, the number of fringes lines passing over a given point {x, y) can 

be calculated for 9 increasing or decreasing from zero to a given wave front angle. 

= y (1 ~ cos 9) — xs\n9 — 2Ad 
^ ^ (A28a) 

y(l-cos0) xsm9 2Ad 

= x A Y (A29) 

The y term in the above equation is positive for 0 ^ 0 and is symmetrical in shape as a function of Q. 
For small angles sin 0 — 0, therefore the coefficient of x is a linear function of angle 0. 
The polarity of rif is therefore dependent on the magnitude and sign of x, 9 and Ad. 
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When 0 — 0, Equation (A29) reduces to 

-2Ad = UfX (A29a) 

Uf is the number of fringes counted at a point (x,y) as Q is varied. Assume Ad = 0, there is a 
special case in Equation (A29) when: 

y(l — COS0) xsxxvd 



X 



(A30) 



and the fringe count rif = 0 despite 0 > 0. In this case, fringe lines would have moved over point 
ix, y) in one direction as 9 is increased and then back again as 9 is increased further to the angle 9 that 
satisfies Equation (ABO). 

As discussed with Equation (A28), the slope of the fringe lines is given by: 



y sin 9 (6 
X (1 — cos 9) 



= COt(j) 



(A31) 



where - is the angle the normal of mirror M2 relative to the j-axis (Figure 3). Therefore, there is a 
set of points {x, y) coincident with the mirror normal that renders Uf = Q. 

Deriving the speed of fringe movement 

If Ad is considered static (A J = 0), Equation (A29) reduces to 

y(l — COS0) xs\n9 
nf = J ^ (A32) 

Taking the derivative of the x-term delivers the speed of concertinaing of the fringe lines at the 
point (x,y). 

^""f - (A33) 



d sin 0 A 

The speed of sideways deflection of the fringe lines at point (x, y) results from calculating the 
derivative of the y-term 

^ =^tane' (A34) 



dsin^ AVl -sin2 0 A 

The overall fringe movement speed is therefore 

dUf X y 



d sin 0 A A 

From Equation (A35), the direction of fringe movement reverses if 



= -; + 4tan0 (A35) 



-Y + Ttan0 = O (A36) 
A A 



And therefore the angle of movement reversal Gygy is 

^ ■ -1 
- = sm ^ -= 

y .JWTy 



.X X 

9rev = tan"^ - = sin~^ , (A37) 
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That is, when the perpendicular of the wave front from the origin points towards +x. The mirror tilt 

9 71 

angle is therefore If y ^ 0, O^gy -, i.e., the smaller y, the larger is 6rev 

From Equation (A29), the fringe count returns to zero if the numerical value of x- and _y-terms 
cancel each other out. 



y — y cos 9 — xsinO — 0 

X 

1 sin 6 = cos d 

y 

1 - -sine = Vl - sin2 6 

y 

X^ X 

1 + ^sin2 0 - 2-sine = l-sin^^ 
y^ y 

X?" X 

— sin^ 0 — 2 — sin 0 = — sin^ Q 

y^ y 



X' X 

-^smQ — 2 — = — sin 0 

y^ y 



x'- X 
^sin^ + sine = 2- 

y^ y 



' x^ 



X 



sin e ^ + 1 = 2 - 



sine = 



y 

2xy^ 



2xy 



This occurs at the angle 0„=o 



e„=o = sin ^ 



2xy 



If X = y, then 0^^=^ = ^, if x > y, then 6^=0 > ^ and conversely , if y > x, then e„=o < \- 
The relationship between 0„=o and Srev results in the identity 0„=o = ^^rev 

2xy 



From Equations (A37) and (A40) 

-1 ■ -1 



X 



= sin 



Resulting in 



^^2 + y2 

2xy 



sirr '^^Y'G'u ^ sin ^^gu cos ^^gu /" 

(x'^ + y^ j 



(A38) 
(A38a) 
(A38b) 
(A38c) 

(A38d) 

(A38e) 

(A38f) 

(A38g) 

(A39) 



(A40) 



(A41) 



(A42) 



And proving the identity of 0„=o = ^^rev being correct. 
Determination of the damping function of the flux curve 

From Figure 4 it is evident that the normalised flux curve oscillates about a level of 0.5 flux and 
decreases its wave amplitude about this level with increasing wave front angle 9. As this behaviour 
constitutes a damped function, the flux decay function is derived subsequently. 
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For a centred photodetector of width s and positioned at x = 0, with integral boundaries X2 — — & 
Xi = — ,y-0 and Ad = 0, the radiant flux ^e(xi,x2) Equation (A8) reduces to 



fks sin G\ ks sin 9 
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1 r //cs sin 0\ 
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ksmdl \ 2 J 



oc 



2 sin J 
k sin 6 



(A43) 



+ s 



In Equation (A43), the maximal flux corresponds to 2s (cf. Equations (All) and (A23): 
flux = constant * z * 2^; z is considered unity as the z-direction is perpendicular to fringe lines). 
Normalising the flux to 2s delivers the normalised flux <I)„ 
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(A44) 



As the normalised flux ranges from 0 to +1, it is converted to a range from -1 to +1 in order to 
compare it to a standard sine function 
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(A45) 



For small 9, Equation (A45) becomes 

/ks6\ 

<D,ocsinc(^— j ^^45^^ 

This equation defines a damped sine function. The equivalent undamped sine function has the form of 

A = sin(2nf6) (A45b) 

Pi 

Where A is the amplitude and /is the reciprocal value of 2nd node angle (= 2-), i.e., the wave 



length of the flux function. Substituting this term into Equation (A45b) yields 

/ s \ /nsO\ 
^ = sm(2^-e) = sm(— ) 



(A45c) 
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The decay function Fq is obtained when normalising the damped sine function to its undamped 
counterpart 

. fks sin d\ 
^ sin(^— ) 

P _ sk sin 6 
^ . fns9\ 

^^"(^j (A45d) 
ks sin 



„ . (ks smU\ 
2sm(^— J 

s/c sin 0 sin (~^) 
After simplifying and considering that 6 is small 



. (ks6\ 
2 sin(^^j 

. (Ttse\ (A45e) 



Considering that k = ^ 

A 



^ _ 2 sin (—J 

° ^'^esin(^) (A46) 
_ 2 1 _ A 1 
sk 6 ns 6 

Equation (A46) is the constitutive equation of the flux decay function at ;c = 0 and j = 0. is a 

2 A 

reciprocal function of G, with a decay constant Co of — or — . 
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